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1.   Introduction 

This  paper  pertains  to  a  device  for  generating  useful  elec- 
trical power  by  means  of  an  electrohydrodynamic  (EHD)  process. 
The  device  utilizes  the  flow  of  a  fluid,  normally  a  gas,  such  as 
air  or  water  vapor,  in  which  are  entrained  a  very  large  number  of 
very  fine  and  well  distributed  solid  or  liquid  particles,  for 
example,  water  droplets  as  in  an  aerosol  spray.   The  particles  are 
electrically  charged,  either  positively  or  negatively,  as  may  be 
convenient  in  a  particular  case.   The  gas  is  caused  to  flow  through 
a  nozzle-like  channel  by  the  imposition  of  a  suitable  pressure  drop. 
Let  station  0  be  the  inlet  to  the  nozzle.   The  charged  particles 
are  introduced  into  the  stream  of  gas  by  a  suitable  injector  at 
some  upstream  location,  call  it  station  i,  and  are  removed  from  the 
stream  by  a  suitable  collector  grid  at  some  downstream  location, 
call  it  station  2.   An  important  intermediate  location,  station  1, 
is  at  the  throat  of  the  nozzle.   The  electrical  charges  which  are 
transported  downstream  along  with  the  particles  on  which  they  re- 
side constitute  an  electrical  current.   These  charges  move  through 
an  electrical  field  which  exerts  forces  upon  them  in  a  direction 
and  sense  opposed  to  the  general  fluid  motion.   By  proper  design  of 
the  injector,  the  size  of  the  water  droplets  or  particles  may  be  so 
regulated  that  these  droplets  have  low  mobility,  that  is,  they  move 
approximately  with  the  surrounding  gas  with  negligible  relative  slip 
produced  by  the  opposing  electrical  forces. 

The  gas  stream  does  work  on  the  charged  particles  in  moving 
them  against  the  resistance  of  the  electrical  forces.   In  this 


process  the  gas  stream  undergoes  a  corresponding  decrease  in 
total  enthalpy.   The  work  done  upon  the  charged  particles  creates 
a  difference  in  electrical  potential  between  stations  i  and  2. 
These  stations  are  connected  by  an  external  circuit  which  includes 
a  useful  electrical  load.   Thus  the  enthalpy  drop  of  the  gas  is 
ultimately  converted  into  a  useful  electrical  power  output  from 
the  external  circuit. 

The  above  scheme  accomplishes  the  primairy  conversion  of  ther- 
modynamic energy  directly  into  electrical  power  without  the  use 
of  any  major  rotating  or  reciprocating  components  such  as  large 
turbines  or  electrical  generators.   Nevertheless,  there  still 
exists  a  requirement  for  a  small  pump  to  recirculate  the  condensate 
and  a  fan  to  recirculate  the  carrier  fluid.   A  condenser  and  a 
boiler  are  also  normally  a  part  of  the  complete  system,  but  these 
various  auxiliaries  are  not  considered  here  in  any  detail  as  our 
present  concern  is  primarily  with  the  nozzle  in  which  the  basic 
thermo-electric  power  conversion  occurs. 

The  above  concepts  are  well  known  and  can  be  found  in  the 
technical  literature.   See  References  1,  2  and  3,  for  example. 
They  are  reviewed  in  this  introduction  merely  to  provide  a  proper 
background  for  our  result  which  is  quite  specific  and  which  is  a 
consequence  of  the  following  analysis. 

The  performance  that  can  be  achieved  by  an  EHD  device  is 
limited  by,  among  other  things,  the  maximum  electrical  field 
strength  that  can  be  sustained  at  the  most  critical  point  in  the 
field  without  inducing  electrical  breakdown,  through  a  spark 
discharge,  of  the  carrier  medium.   The  electrical  field  that  is 


present  in  most  EHD  devices  has  its  highest  value  at  the  point 
of  injection.   See  References  1-4,  for  example.   In  this  regard 
there  is  experimental  evidence  that  breakdown  field  strength  in 
any  region  varies  linearly  with  the  quantity  Rp  where  R  is 
the  gas  constant  and  where   p   is  the  fluid  density  in  the  region 
considered.   This  experimentally  established  limit  depends  of 
course  on  the  type  of  gas/aerosol  combination  that  is  involved. 

But  the  presence  of  a  maximum  field  at  the  beginning  of  the 
EHD  conversion  channel  limits  performance  severely  because  the 
rest  of  the  channel  must  operate  below  maximum  capacity  or,  in 
other  words,  the  electric  pressure  is  at  its  allowable  limit  in 
only  a  small  fraction  of  the  conversion  channel. 

The  above  considerations  suggest  that  an  optimum  design 
would  be  one  in  which  the  local  field  strength  is  everywhere  uni- 
formly close  to  the  critical  limit.   Such  a  design  would  achieve 
greater  electric  work  output  per  unit  mass  of  fluid  than  would 
any  other.   The  essence  of  this  paper  is  the  recognition  that  a 
uniform  maximum  loading  through  the  channel  is  the  optimum  load- 
ing and  that  this  condition  can  indeed  be  achieved  at  least  to 
an  acceptable  approximation.   A  secondary  aim  of  this  paper  is - 
to  derive  the  characteristic  geometrical  features  of  the  optimum 
design  and  certain  associated  performance  parameters  and  limits. 

The  geometrical,  electrical,  and  thermodynamic  features  of 
such  a  power  generator  are  governed  by  various  physical  laws  of 
which  one  of  the  most  significant  in  the  present  context  is 
Poisson's  equation  as  it  applies  to  an  electric  field.   In  its 
full  generality,  Poisson's  equation  is  three  dimensional  in 


nature;  by  restricting  attention  to  configurations  having  polar 
symmetry  we  can  simplify  this  to  a  two-dimensional  form.   More- 
over, by  further  restricting  the  application  to  configurations 
whose  largest  radial  dimensions  are  small  in  comparisons  with 
their  axial  length,  we  may  finally  simplify  Poisson's  equation  to 
a  one-dimensional  approximation.   Despite  certain  limitations  in 
accuracy,  this  one-dimensional  version  can  be  more  enlightening 
than  the  more  elaborate  two-dimensional  analysis  because  it  shows 
basic  trends  so  much  more  simply  and  clearly.   It  is  this  one- 
dimensional  approximation  on  which  the  present  paper  is  based. 

The  analysis  must  deal  not  only  with  Poisson's  equation  for 
electrical  fields,  but  also  with  the  laws  of  fluid  flow.   It  is 
consistent  with  the  foregoing  commentary  to  write  these  laws  also 
in  their  common  one-dimensional  forms.   At  this  level  of  ideal- 
ization it  is  also  appropriate  to  treat  the  flow  through  the 
channel  as  isentropic.   We  also  consider  the  fluid  as  a  perfect 
gas  for  which   y  /  the  ratio  of  specific  heats,  is  constant. 

One  other  special  assumption  is  involved.   This  is  based 
on  the  fact  that  the  drop  in  total  temperature  through  the 
channel  is  normally  very  small  in  comparison  with  the  absolute 
total  temperature   T   at  the  nozzle  inlet,  station  0.   Hence  it 
is' permissible  for  certain  analytical  purposes  to  neglect  this 
temperature  drop  and  to  approximate  the  flow  through  the  channel 
as  an  isentropic  flow  of  constant  total  temperature.   Calculations 
to  be  presented  later  show  that  the  actual  change  of  total  tem- 
perature amounts  to  less  than  one  percent  for  a  typical  case. 


Future  progress  in  EHD  power  generation  might  increase 
electric  power  output  by  an  order  of  magnitude  over  levels  that 
seem  currently  feasible.   In  that  case  the  drop  in  total  tempe- 
rature through  an  EHD  nozzle  could  become  significant.   Fortunate- 
ly, the  present  analysis  can  be  revised  and  generalized  to  account 
rigorously  for  such  variation  in  total  temperature  when  and  if 
necessary.   Under  present  circumstances,  however,  it  is  not 
warranted  to  complicate  the  analysis  to  include  this  refinement 
since  at  low  power  outputs  it  has  no  appreciable  effect  on  the 
final  calculated  results. 
2.   Analysis 

It  is  convenient  for  the  purposes  of  the  present  analysis  to 
introduce  a  parameter   a  which  is  used  to  distingiiish  between 
the  two  distinct  cases  of  positively  and  negatively  charged  par- 
ticles.  Specifically,  we  set  a  =  +  1   for  the  case  of  positively 
charged  particles,  and  we  set  a  =  -  1   for  the  case  of  negatively 
charged  particles.   It  is  then  appropriate  to  denote  the  electrical 
charge  per  unit  mass  by  means  of  the  product  aq  where   q   has 
units  of  coulombs/kg  and  is  always  positive  by  definition  while 
a   is  a  dimensionless  factor  as  defined  previously. 

Using  the  above  notation  we  may  write  the  following  two 
expressions  for  the  electric  current  i   and  the  mass  flow  rate 
m   through  the  channel,  namely, 

i  =  aqpAv  =  constant  along  channel  (2.1) 

ih  =  pAv   =  constant  along  channel  (2.2) 


where 

i  =  electric  current,  amps 

ih  =  mass    flow   rate    through   channel,    kg/sec 

q   =  electric   charge   per   unit  mass,    positive  by 

definition,    coulomb/kg 
p    =   fluid  density,    kg/m      (variable) 

A  =   cross-sectional    area   of   channel   at   any 

2 
given   strearawise   station,      m      (variable) 

V  =  mean  axial  velocity  at  any  given  streamwise 
station,  m/sec  (variable) 

Note  that  Eq.  (2.1)  is  based  on  the  previously  stated 
assumption  that  size  of  the  fluid  particles  and  the  charge  on 
each  particle  can  be  so  regulated  that  the  mobility  of  the  par- 
ticle with  respect  to  the  surrounding  medium  is  essentially 
negligible. 

By  dividing  Eq.  (2.1)  by  Eq.  (2.2)  we  also  find  that 

-r-  =  aq  =  constant  along  channel  (2.3) 

m    ^  ^ 

The   one-dimensional    version   of  Poisson's    equation    can    now 
be  written   in    the    form 


(2.4) 


where   z  =  axial  coordinate,   m 

<t)   =  electrical  potential,  volts 
£  =  £   =  electrical  permittivity  of  free  space 
=  8.854  X  lO"^  farad/m 


The  local  electrical  field  strength  at  an  arbitray  axial 
station  z  is  defined  as  usual 

^  =  -  (§)  '2-5' 

where  E  =  field  strength,  volts /m 

As  will  be  seen  here,  an  optimum  EHD  generator  has  an  elec- 
tric field  which  does  not  change  sign  from  inlet  to  outlet. 
Under  these  circumstances  the  electric  field  is  always  negative 
for  positive  space  charged  and  positive  for  negative  space  charge. 
Hence  we  may  write 

E  =  -  a|E|  (2.6) 

It  has  been  established  by  experiment  that  over  a  broad  (but 
not  unlimited)  range,  the  field  strength  at  breakdown  is  well 
represented  by  the  simple  linear  law  (see  Reference  2) . 

I^bI  =  =o  -^  S  '^  P  '2.7) 

where   C   and  C_   are  characteristic  constants  of  the  medium, 
o        B 

R  is  the  gas  constant  in  units  of  Joule/kg  K   .   Experimental 

measurements  show  that  constants   C   and  C„   happen  to  have  the 

o        B    '^^ 

same  numerical  values  for  both  air  and  steam.   However,  we  defer 
to  later  pages  of  this  paper  any  reference  to  actual  experimental 
values.   Of  course,  the  present  analysis  is  restricted  to  that 
specific  range  of  conditions  for  which  Eq.  (2.7)  is  in  fact  a 
valid  approximation. 

If  we  now  impose  the  constraint  that  the  local  field  strength 
is  everywhere  just  equal  to  its  critical  value  at  impending  break- 
down, we  readily  find  from  the  last  three  equations  that 


II  =  alC^   +  =8  ""'  '2.8) 


Now  differentiating  this  gives 


=  a   r      R  ^  (2.9) 


Upon  equating  the  right  sides  of  Eqs.  (2.4)  and  (2.9)  and 
simplifying,  we  obtain  the  important  relation 

i  (^1  = 2_  (2  10) 

The  reciprocal  of  the  quantity  on  the  right  side  of  Eq.  (2.10) 
(momentarily  disregarding  the  negative  sign)  now  identifies  a 
significant  characteristic  length;  let  us  denote  it  by  symbol  X    . 

£  C   R 
X  = —    (meters)  (2.11) 

q 

It  is  very  instructive  to  rewrite  Eq.  (2.11)  in  the  follow- 
ing alternative  form 

Xq  =  X  .   q„^^  =  £  C_  R  =  constant  (2.12) 

^    mm   max      B 

Notice  that  the  quantity   £  C   R  is  a  characteristic  pro- 
perty  of  the  medium.   Once  the  medium  is  chosen  the  designer  has 
no  further  control  of  the  value  of  this  constant.   Eq.  (2.12) 
now  tells  us  that  in  order  to  maintain  the  electrical  loading  at 
incipient  breakdown,  the  product   Xq  must  remain  constant.   Thus 
any  increase  in  q   must  be  accompanied  by  a  corresponding  decrease 
in   X  ,  or  vice  versa. 


Experience  up  to  now  indicates,  however,  that  there  exists 
some  practical  upper  limit  q^i^^   ^°  ^®  value  of  q   that  can 
actually  be  achieved  for  any  specific  type  of  design,  and  hence 

there  exists  some  corresponding  lower  limit   X_ .    on  the  accom- 

^  ^  mm 

panying  longitudinal  characteristic  length.   Of  course,  a  value 

of  q   less  than  q     can  always  be  employed  if  necessary  in 

which  case  the  corresponding  value  of   A  will  be  greater  than 

X  .   . 
mm 

The  situation  is  shown  schematically  in  Fig.  2.1  in  the  form 
of  a  log-log  plot.   The  solid  line  represents  the  locus  of  points 
all  of  which  produce  incipient  breakdown.   Incidentally,  it  will 
also  be  shown  later  that  all  of  these  points  correspond  to  the 
same  theoretical  power  output  from  the  generator.   Note  that  the 
line  continues  on  indefinitely  toward  decreasing  values  of  q 
and  increasing  values  of   X  .   At  the  right  end,  however,  it  ter- 
minates at  the  point  corresponding  to  q     and   X  .  . 

^  r     i.     ^max        mm 

Some  typical  values  of  these  various  characteristic  constants 
based  on  p\±ilished  data  are  summarized  in  Table  2.1  for  reference 
purposes . 

Upon  substituting   X   into  Eq.  (2.10),  we  can  easily  inte- 
grate the  simple  expression  that  is  thereby  obtained.   The  result 
is 


(2.13) 


Locus  of  Possible  Design  Points  for 
Critical  Electrical  Loading 


'max 


q    (log  scale) 


Fig.    2.1      Design   Points    and   Design    Limits 

for   an   OptimiJiin   Slender   EHD  Nozzle 
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Table    2.1      Typical  EHD  Properties 


Source  o 

f  Data 

kef.  1 

Ref.  2 

Ref.  3 

Medium 

Water 
droplets 
in  steam 

Water 
droplets 
in  air 

Water 
droplets 
in  air 

Quantity 

Units 

^o 

volts/m 

8.63  X  10^ 

8.63  X  10^ 

8.63  X  10^ 

S 

m^  Vc 

9.49  X  10^ 

9.49  X  10^ 

9.49  X  10^ 

e  Cg  R  =  Aq 

C  m/kg 

3.9  X  10"^ 

2.4  X  lO"^ 

2.4  X  10~^ 

^max 

C/kg 

2.5  X  10"^ 

3  X  10"^ 

3  X  10"^ 

X  . 
min 

m 

1.5  X  10"^ 

8  X  lO"^ 

8  X  10""^ 

This    important  result  is    the   key    to   all   remaining   details   of 
the  analysis. 

Eq.     (2.13)    reveals    that   the    density   continues    to    drop  mono- 
tonically  with   increasing   distance   downstream.      Thus    the   flow  cor- 
responds   to   the   known    isentropic    flow   through   a    converging-diverging 
nozzle    for  which    the  walls    are    so    designed   that    the   density   drops 
exponentially   as    specified  by  Eq.     (2.13) .      Let  subscript    1   denote 
the   throat  of   the  nozzle.      Then   the    following   relations    are  known 
to   apply 


l.^L^M^ 


(y  -   1) 


(2.14) 
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r 
r. 


1  + 


r-^M^^ 


Y  +  1 


(Y 

4- 

1) 

4(y 

1) 

1 

(2.15) 


where   y     =  ratio  of  specific  heats 
M  =  Mach  number  at  station  z 
r  =  nozzles-radius  at  station  z 


r,  =  nozzle  radius  at  throat 


It  is  readily  deduced  from  Eqs.  (2.13)  and  2.14)  that 


(f)  = 


(Y  -  1) 


In 


l.I-^M^ 


(2.16) 


Eqs.  (2.15)  and  (2.16)  now  constitute  a  pair  of  parametric 

r  z 

equations  for   ( — )   as  a  f-onction  of  i^)      with  Mach  number  M 

r^  A 

as   parameter.      These   equations    define    the    shape   of    the   optimum 

nozzle   having  uniform  loading   just  below  breakdown   all   along  its 

r  z 

length.      A   tabulation   and  a   curve   of      ( — )       versus       (-r-)       are 

r^  A 

given  in  the  next  section,  using   A  =  1.3  .   This  value  of   A   is 
representative  for  steam.   See  Table  4.1  and  Fig.  4.1. 

At  the  throat  of  the  nozzle   M  =  M  =  1  .   Then  Eq.  (2.16) 
gives 


I 


(Y  -  i: 


£n 


Y  +  1 
2 


(2.17) 


For  y   =   1.2      this  yields 


-^  \   =  0.466 


(2.18) 
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Let   z.   denote  the  axial  station  which  locates  the  injector 
that  introduces  unipolar  electrical  charges  into  the  flow.   We 
do  not  stipulate  in  advance  any  particular  fixed  location  for  the 
injector  because  the  nature  of  the  charged  particles  may  require 
one  location  or  another  depending  on  certain  supersaturation  con- 
ditions that  might  be  needed  for  the  formation  of  the  aerosol 
droplets  (see  Reference  1) .   Moreover,  the  converging-diverging 
nozzle  does  not  have  to  be  of  the  shape  specified  in  this  analysis 
upstream  of  the  injector  location. 

Of  great  importance  is  the  theoretical  gross  power  output  of 
such  an  optimum  device.   To  determine  this  we  must  first  find  the 
current  from  Eq.  (2.1).   It  is  advantageous  to  evaluate  the  terms 
on  the  right  side  of  Eq.  (2.1)  at  the  nozzle  throat,  station  1. 
Thus, 


i  =  aqp^A^v^ 


(2.19) 


Now  letting   a   and  a.   denote  sonic  velocity  at  stations 
0  and  1,  respectively,  we  have 


v^  =  a,  =  a 

1    1    o 


(2.20) 


Upon  combining  the  above  two  relations  we  finally  obtain  the 
current  in  the  form 


i  =  aqp  a  A^ 
^  o  o  1 


(2.21) 
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Note,  however,  that  from  the  previous  analysis  of  the  isen- 
tropic  nozzle  flow,  it  follows  that 


Pl\Pl 


Y  +  1) 

n  2(Y  -  1) 


Y  +  1 


(2.22) 


Moreover,  it  is  useful  to  eliminate  q   from  Eq.  (2.21)  by 

means  of  Eq.  (2.12)  and  to  rearrange  the  result  in  dimensionless 

form.   In  this  way  we  finally  obtain  the  dimensionless  current  in 
the  form 


Xi 


£  C_  R  p  a  A- 
B    ^o  o  1 


=  a 


Y  +  1 


(Y  +  1) 

2(y  -  i) 


(2.23) 


1 


The  next  step  is  to  find  the  change  in  electric  potential 
between  stations  i  and  2.  This  is  accomplished  by  integrating 
Eq.  (2.8)  between  these  limits.  We  again  rearrange  the  result 
in  dimensionless  form.   The  result  is 


z . 

1 


(^2  -  *i) 
Po  S  ^  ^ 


=  a 


"2  -  h 


,S  ""   '. 


-  e 


(2.24) 


Finally,  the  theoretical  power  output  follows  from  the  rela- 


tion 


P^  =  i($2  -  *i) 
where   P   =  theoretical  gross  electrical  power,  watts 


(2.25) 
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The  desired  power  can  be  obtained  in  dimensionless  form  by 
direct  multiplication  of  Eqs .  (2.23)  and  (2.24).   After  some  minor 
rearrangement  we  obtain 


(y  +  B 

2(Y  -  1) 


^  2  „  2   2    T 

£  Cg  R   p^  a_  A. 


o  IJ 


-  e 


Y  +  1 


(2.26) 


This  result  is  fundamental.   It  shows  just  how  the  theore- 
tical gross  power  output  of  the  ideal  nozzle  depends  on  the  various 
geometrical  and  physical  parameters  of  the  problem. 

It  is  also  instructive  to  observe  from  Eqs.  (2.23),  (2.24) 
and  (2.26)  how  current,  voltage  and  power  vary  with  the  character- 
istic length  parameter  X      of  the  nozzle.   In  other  words  we  hold 
nozzle  shape  constant  and  simply  vary  the  absolute  size.   Thus 
the  quantities   z./X   and  "^2^^      remain  constant.   Then  over  the 

range  X   >    X    .       , 
^  —     mm  ' 


i  ~  X 


-1 


(2.27) 


(^2  -  \^ 


(2.28) 


P  ~  X   =  independent  of  X 


(2.29) 


Hence  current  is  inversely  proportional  to  size,  voltage 
is  directly  proportional  to  size,  and  power  is  independent  of 
size'.   This  comes  about,  of  course,  because  of  the  limitations 
imposed  by  the  electrical  breakdown  phenomenon'  see  also  Ref.  2. 
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Voltages  tend  to  be  inconveniently  high  in  EHD  devices . 
The  above  scaling  rules  suggest  that  this  problem  can  be  alle- 
viated by  making  EHD  devices  small,  with  X      as  close  to  X    . 

mm 

as    feasible.      To    the   extent   that  it   can  be   practically  achieved, 
this   method  has    the   advantage    that  it   decreases    the   size   of   the 
device  without  any   corresponding  decrease  in  its   power  output. 

Incidentally,  it  should  not  be  overlooked  that  Eq .  (2.26), 
being  based  on  a  purely  one-dimensional  approximation,  has  high 
accuracy  only  for  sufficiently  low  values  of  the  ratio  ^9^ ^2  * 
3.      Optimization 

It  is  convenient  to  simplify  Eq .  (2.26)  by  introducing  the 
following   auxiliary   notation.      Let 


?=x  "-^i 


and   let  Eq.     (2.26)    be    rewritten   in   the    form 


P 


^  =  F   = 


.^  S'  ^  '  Po'  ^o  h^        ' 


-C.  -  Col 


a(^2  -  ?i)   +  e       ^  -  e 


fy  +  D 
2(T  -  1) 


Another  quantity   of    fundamental    importance   is    the  exit 
kinetic  power  of    the    jet.      This    is    defined   as    follows 


V2^ 


By  using  the  various  well  known  relations  for  isentropic 


16 


compressible  flow  in  a  channel/  we  can  readily  reduce  this  to 
the  following  dimensionless  form 


P  . 


P   a  ^  A, 
^o   o    1 


=  F.CM^) 


(3.5) 


where  function   F . (M_)   is  as  follows 

J      2 


For     M      <    1 


(3y   -    1) 
2(Y    -    1) 


M   ■ 
^j^^2'  2 


1    .      S^      M,^ 


(3.6) 


For     M-   ^  1 


F.  (M.) 

:      2 


1 

2 

(Y    +    1) 

2(y  -  1) 

2 

Y    +    1 

-  1 


1  + 


Y    -    1 


M, 


Mj        (3.7) 


At     M2   =5^  1    -these    two   expressions   both   give    the   result 


(3y    -    1) 

2(Y    -    1) 


F.(l) 


2 

Y 

+ 

1 

(3.8) 


It  will    also  be   useful   later   to  have    the   value   of      F. 

1 

in  the  limit  as   M2   increases  without  bound.   From  Eq .  (3.7) 
we  may  deduce  that  this  is 


(Y  +  1) 

2(y  -  1) 


Fj(«) 


(Y  -  1) 


Y  +  1 


(3.9) 
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Consider  an  EHD  power  generator  which  operates  as  a  closed 
system.   The  fluid  which  discharges  from  the  nozzle  exit  must 
be  recirculated  back  to  the  nozzle  inlet.   This  involves  dece- 
lerating the  flow,  cooling  the  fluid,  condensing  and  separating 
the  liquid  droplets,  and  so  on.   It  involves  restoring  the 
pressure  of  the  fluid  by  means  of  a  suitable  fan,  compressor  and/ 
or  pump,  reheating  the  liquid  and  completing  the  thermodynamic. 
cycle.   These  processes  require  power.   The  basic  compressor  and/ 
or  pump  power  requirement  is  inherent  in  the  thermodynamic  cycle 
itself  and  cannot  be  eliminated.   Additional  power  demands  arise 
from  the  various  irreversibilities  that  inevitably  occur  in  any 
real  physical  system;  this  additional  power  cannot  be  completely 
eliminated  but  it  can  be  minimized  by  careful  engineering  design. 
Fortunately,  some  of  the  power  required  can  actually  be  recovered 
from  the  exit  kinetic  power  of  the  jet  itself.   Inevitably,  how- 
ever, there  will  be  a  net  power  demand,  call  it   5  P .  ,  which  must 
be  subtracted  from  the  gross  electric  power  output  P   .   Conse- 
quently, the  net  useful  electrical  power  becomes 

P   ^^^  =  P   -  5  P.  (3.10) 

e  net    e      j 

Yet  another  jet  power  term  is  very  useful.   This  is  the  jet 
power  developed  at  the  throat  of  the  nozzle,  call  it  P.(l)  . 
This  makes  a  useful  dimensional  reference  parameter  which  may  be 
used  to  normalize  the  net  power  output  P       . 

This  reasoning  leads  us  to  define  the  following  useful  net 
power  coefficient,  namely. 
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P        .  (P^  -  <s  p.) 

e   ^et  =   Y  =        Q  1-  (3   11) 

P.  (1)         ^  P.(l)  VJ.i-i-; 


Upon   substituting   into   this   definition   the   expressions    for 

P      ,    P.      and     P.(l)      as   previously   developed,    and  upon   simpli- 
®  J  J 

fying,   we   can   reduce   the   result  to   the   form 


(3.12) 


where,    for  a  perfect  gas, 

It   is    therefore    clear   that  we   should  maximize    inlet  stagna- 
tion pressure      p        while   keeping      T        no  higher    than   necessary    to 
obtain  proper   droplet  size   in   the  medium.      Moreover   it  might  be 
possible    to   eliminate   this    restriction  on      T        by   injecting   solid 
or   liquid   charged  particles    rather    than    condensing  droplets.      In 
this    case      T        could  be   sharply   reduced   thus    greatly   incresing 
the   value   of   the   crucial  parameter      3    . 

Eq.     (3.12)       for   the   power   coefficient      x      is    the   key   result 
of   this    section.      The   optimization   of    an  EHD   generator   can  be 
based  in  part  on   an  effort   to  maximize    this   quantity.      In    this 
connection,    however,    another  quantity  of   interest  is    a  power 
ratio      r\      which  we   define   as    follows : 

P  .\  5    F. 

^^  Un  =  l-gT:i  (3.14) 

e       /  ^      e 
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Notice  that  the  quantity  ri   as  defined  above  does  not 
allude  to  the  net  heat  input  to  the  cycle.   Hence  it  should  not 
be   confused  with  overall  thermal  efficiency  which  will  always 
be  some  very  much  lower  figure. 

The  various  power  terms  that  occur  in  Eq.  (3.12)  can  also 
be  used  to  determine  the  ratio  of  AT   ,  the  drop  in  total  tem- 
perature through  the  nozzle,   T   ,  the  inlet  total  temperature. 
This  works  out  to  be  simply 

AT     6  F 

9.   =     g  (3  15) 

T      F.  (0°)  l-^i.lii 

O      3 

Recall  that  the  present  analysis  is  restricted  to  cases  in 
which   AT  /T   is  small  compared  with  unity.   Eq.  (3.15)  may  be 
used  to  check  whether  this  restriction  is  in  fact  satisfied  in 
amy  particular  case. 

It  is  evident  from  Eq.  (3.12)  that  the  value  of  the  net 
power  coefficient   x   that  can  be  attained  is  governed  by  the 
two  key  parameters   S   and   6  .   The  coefficient   6   controls 
the  maximum  gross  power  that  can  be  generated  without  electrical 
breakdown,  and  the  factors  that  govern   S   in  turn  are  shown  in 
Eq.  (3.13) .   Coefficient   6   describes  the  power  loss  inherent 
in  the  thermodynamic  cycle  itself,  as  augmented  by  the  various 
irreversibilities  that  must  inevitably  occur  in  such  a  system. 
4.   Summary  of  Key  Relations  and  Results 

The  maximum  specific  charge   q   that  can  be  employed  in  an 
ideal,  slender  EHD  nozzle  without  causing  electrical  breakdown 
is  governed  by  the  relation 


20 


Aq  =  X  .   q    =  e  r  R  =  constant 
^    min  ^max      B 


(4.1) 


The  proper  geometric  form  for  an  optimum  EHD  nozzle  is 
defined  by  the  following  pair  of  parametric  equations 

(y  +  1) 


1  / 


1  +  1^  m2 
2 


^  4(Y  -  1) 


1 

/M 


(4.2) 


^=«X 


(y  -  i) 


^n 


1.^4^m2 


(4.3) 


The  value  of  z,^      corresponding  to  any  chosen  value  of  exit 
Mach  number  M   may  be  calculated  from  Eq .  (4.3) .   Then  we  set 


a  = 


S^^o 


(4.4) 


and  find 


F   = 


a(<;^  - 


=i' 


+  e 


-  ^i 


-  ?. 


-  e 


(Y  +  1) 
2(Y  -  1) 


Y  +  1 


(4.5) 


where   ^.   denotes  the  location  of  the  injector. 

The  following  functions  and  constants  may  next  be  found. 
For  M  <_  1 


F.  =  -^ 

J     2 


From  M_  >_  1 


3  r 


I.iX^Im,^ 


(3y  -  1) 
2(y  -  1) 


(4.6) 
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F.  =  -^ 


2 

(Y    +    1) 
■2(y    -    1) 

Y    +    1 

^  -  1 


1  + 


Y  -  1 


M, 


also 


F.(l) 


(Y    -    1) 

2 

2(Y    -    1) 

[y  +  ij 

and 


F  .  (") 


Finally 


1 
(Y  -  1) 


Y  +  1 


(Y  +  1) 
2(Y  -  1) 


(4.7) 


(4.8) 


(4.9) 


e  = 


'   S  ^o 

Y  T  ^ 


(4.10) 


^  ~  F-TU 


S  F   -  5  F. 


(4.11) 


5  F. 


(4.12) 


AT 


B  F 


F.(«') 


(4.13) 


Values  of   X  /  n   and  AT  /T   versus   M^   are  listed  in 
Table  4.2  for  a  typical  case  which  is  defined  by  the  specified 

The  value   6  =  0.1   corres- 
ponds to  air  or  steam  at  about  250  atmospheres;  this  example 


values  of  y    ,    c, .     ,    a    ,    ^      and   6 
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illustrates  the  extreme  measures  that  are  required  to  lift  per- 
formance to  levels  which  are  even  marginally  acceptable.   The 
assumed  value   6  =  0.1   implies  that  pump  work  and  all  fluid 
losses  amount  to  only  10  percent  of  the  kinetic  power  of  the  jet; 
this  is  a  very  optimistic  hypothetical  example,  but  it  suffices 
to  show  certain  trends.   Reference  (5)  gives  more  detailed  data 
on  one  particular  type  of  loss  that  may  be  involved,  namely, 
losses  in  pipe  or  channel  bends. 

Table  4.2  shows  clearly  the  general  effects  of  varying  exit 
Mach  number  ^^  •   It  also  shows  the  trade-off  between  perfor- 
mance parameters   x  ^^'^     n  '  generally  speaking  an  increase  in 
either  of  these  can  only  be  obtained  at  the  expense  of  a  decrease 
in  the  other.   Note  that  excessively  high  values  of  M   lead  to 
negative  values  of  outputs 

The  last  column  of  Table  4.2  reveals  that,  over  the  range 
of  greatest  interest,  the  value  of  AT  /T   is  smaller  than  0.01. 
This  confirms  that  one  of  the  basic  restrictions  of  the  theoretical 
analysis  is  indeed  satisfied  in  this  case. 

The  optimum  EHD  device  disclosed  herein  is  useful  not  only 
as  a  possible  real  physical  system  but  also  as  a  hypothetical 
standard  against  which  the  performance  of  any  other  real  or  hypo- 
thetical EHD  system  may  be  compared. 

The  analysis  of  this  section  could  be  greatly  improved  by 
treating  the  optimum  EHD  channel  as  just  one  element  of  a  com- 
plete thermodynamic  cycle.   By  working  with  such  a  complete  cycle 
in  detail  the  true  value  of  parameter   6  ,  instead  of  being 
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TABLE   4.1 
OPTIMUM  SLENDER  EHD  NOZZLE 


Y    =    1.3 

M 

z/X 

'/""i 

) .  2  0 

0.0100 

1.  7  3-"'3 

I'.kl 

0.0  701 

1. 25 5 S 

) .  3  0 

0.1753 

1 .  0  0  2  f; 

0  .  5  0 

0  .  '3  0  5  3 

1.013  '.3 

i ,  0  0 

0 .  U  3  5  0 

1. 0  00  0 

1.20 

0.3  310 

1.0150 

l.i+0 

0.3501 

1. 0  5  96 

1.6  0 

1.3833 

1.  1275 

1.33 

1.  32  0  3 

1. 2132 

2.  0  0 

i  .  3  3  a  7 

1. 5315 

2.2  0 

1  .  0  1  J  4 

1.4522 

2.4  0 

2.0737 

1.5200 

2.30 

2.3337 

1. 3153 

2.3  0 

2. 5016 

2.0230 

3.00 

2.  3i+3  1 

2, 2715 

3.20 

3.10  20 

2.5451 

3.4  0 

3.3525 

2.2518 

3.30 

3.5  H2 

3.1338 

3.8  0 

3.8416 

3.573:+ 

i+.  0  0 

'4.  0  70  3 

3.0330 

4.  20 

:+.3121 

k  ,'lr   5   51 

1+ .  '■+  0 

'+.5^00 

1+ .  0  5  2  3 

U  .  '■>  0 

!+.  75  2  0 

5. 5  155 

4.3  0 

4.1:02 

3. 1214 

5.0  0 

5.1038 
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TAB-LE    4.2       PERFORMANCE    OF    OPTIMUM   EHD    GENERATOR 


Y    =   1.3  ^.    =0  a=0  6=   0.1  5    =   0.1 


^2 

4 

X 

ri 

AT  /T 

0.>05 

0*0012 

0*0003 

0.9146 

0  *  0000 

0  ♦  1 0 

0*0050 

0*0010 

0.3297 

0*0001 

0  ,  1  s 

0,0.1.12 

0,0019 

0*7459 

0*0003 

0,20 

0*0199 

0*0030 

0  *  6636 

0.0G0(;> 

0,25 

0,0311 

0*0041 

0  *  5333 

0*0009 

0  ♦  30 

0*0447 

0*0051 

0*5054 

0*0013 

0.35 

0,0607 

0*0058 

0*4304 

0*0013 

0 .  40 

0*0791 

0*0063 

0,3535 

0*0023 

0*45 

0,0997 

0  *  0063 

0*2903 

0^0023 

0,50 

0,1227 

0.0060 

0*2260 

0  *0035 

0*55 

0,1479 

0*0052 

0*1653 

0*0041 

0,60 

0,1753 

0*0041 

0* 1099 

0  *  0043 

0*65 

0*2043 

0*0025 

0>0537 

0^0056 

0*70 

0*2364 

0*0006 

0  *  0 1 2 1 

0*00^3 

0*75 

0*2700 

-0*0016 

-0*0297 

0*0071 

0,80 

0*3056 

-0*0040 

-0*0663 

0*0079 

0 ,  35 

0*  3430 

-0*0066 

-0*0990 

0.0087 

0 .  90 

0*3822 

-0*0092 

-0*1265 

0,0095 

0*95 

0*4232 

-0*0113 

-0*1493 

0*0104 

J.  >  00 

0*4659 

-0*0143 

-0* 1675 

0*0112 
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merely  assiimed,  could  be  actually  calculate  for  each  case.   Also 
instead  of  working  with  relative  indices  of  performance  like   x 

and  n  ,  we  could  calculate  actual  overall  thermodynamic  effi- 
ciency.  Moreover,  the  effect  on  overall  thermodynamic  efficiency 
of  varying  various  characteristics  of  the  thermodynamic  cycle 
could  then  be  systematically  evaluated.   Time  limitations  have  pre- 
cluded such  a  more  ambitious  and  informative  effort.   The  present 
report  concentrates  instead  primarily  on  the  concept  of  the  opti- 
mum EHD  channel  itself. 
5.   Comparison  of  Optimum  Channel  with  Uniform  Channel 

In  EHD  power  generation,  a  uniform  channel  has  the  obvious 
merit  of  simplicity.   In  this  case  the  gas  density  remains  prac- 
tically constant.   However,  the  electrical  field  is  always  a  maxi- 
mum at  the  entrance  to  the  working  section,  and  gradually  decreases 
thereafter.   Maximum  power  output  is  obtained  if  the  length  be  so 
chosen  that  the  electrical  field  is  zero  at  the  exit  from  the 
working  section. 

It  is  well  known  that  the  electrical  power  output  of  such  a 
uniform  channel  is 

P   =  ^  e  E_^  AV      (watts)  (5.1) 

e    2    B 


where 


E   =  breakdown  strength  of  medium,  volts/m 

2 

A  =  area  of  uniform  channel,  m 

V  =  constant  mean  velocity  in  channel,  m/sec 
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For  the  purpose  of  the  present  section  we  divide   P   by 
the  mass  flow  rate  in  thus  fixing  the  electrical  work  per  unit 

mass  of  fluid  (Joule/kg) .   It  is  also  convenient  to  nondimen- 

2 

sionalize  the  resulting  quantity  by  dividing  through  by   a*  /2 

(Joule/kg)  where   a*   is  the  sonic  velocity  which  corresponds  to 
a  Mach  number  of  unity.   We  also  utilize  Eq.  (2.7)  and  make  use 
of  the  fact  that  the  constant  C   is  negligible  in  comparison 
with  the  quantity   ^^P  •   ^^  addition  we  employ  various  standard 
perfect  gas  relations.   In  this  way  we  reduce  Eq.  (5.1)  to  ob- 
tain a  dimensionless  power  coefficient  C   as  follows: 

^  e 


ft  a*^/2     ^      (pAV)  [(^^)  RtJ 


(5.2) 
2  .-  (7^ 


1  +  ^-^r^   m2 


This  shows  clearly  that  reducing  the  Mach  number  M  increases 
the  electric  work  per  unit  mass  of  fluid.   The  maximum  value  that 
can  be  achieved  in  this  way  is 

Now  let  us  consider  the  optimally  loaded  configuration. 
Eq.  (2.26)  provides  a  convenient  starting  point.   As  in  the  pre- 
vious case  we  neglect  the  quantity   (C  /C„  R  p  )   thus  eliminating 

oso 

the   corresponding    term   from   the  equation.      The  maximum  power   is 
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obtained  by  setting  z./\   =   0   and  z^/X  =  °°    .   The  result  can 
then  be  nondimensionalized  in  much  the  same  way  as  for  the  pre- 
vious case.   Notice  that  in  this  case  sonic  velocity  occurs  at 
the  throat  of  the  nozzle  so  that  a*  =  a.  .   The  required  result 


can  now  be  developed  as  follows 


m  a*  /2 


=  C 


opt 


opt 


/   ^  2  „2    2 

^^  S   ^   PQ   ^o  \' 

(p^  A^  a^)  (a^V2) 


(y  +  1) 
2(y  -  1) 


Y  +  1 


(Y  +  1) 
2(Y  -  1) 


^  S'  ^^Po)^(?f^(^)^ 


C: 


4—1 )  RT 
Y  +  1    .o 


Y  +  1 


(5.4) 


1     ,  1     (Y  +  1) 
(Y  -  1)    2  "    2(Y  -  1) 


-  ^^  i^MU^  ^oi 


Y  +  1 


The  bracketed  exponent  turns  out  to  be  zero  whereupon  the 
final  result  becomes  simply 


m  a*^/2 


=  C 


opt 


opt 


/y  +  1 
Y 


'z    C 


B   ^o 


(5.5) 


Inspection  of  Eqs.  (5.3)  and  (5.5)  now  reveals  that  the  op- 
timally loaded  channel  produces  twice  as  much  electric  work  per 
unit  mass  of  fluid  as  does  the  uniform  channel.   These  formulas  ^ 
are  of  fundamental  importance  because  they  show  the  basic  per- 
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formance  limits  of  these  two  types  of  EHD  channels.  Substitu- 
tion of  typical  numerical  values  into  these  formulas  will  also 
show  how  stringent  these  performance  limits  really  are. 

An  improvement  in  performance  by  a  factor  of  2  is  unfortu- 
nately not  sufficient  in  itself  to  bring  EHD  power  generation  to 
practical  reality.   A  further  improvement  by  a  factor  of  somewhere 
between  5  and  50  is  probably  still  needed  to  achieve  that  goal. 
Moreover  a  factor  of  2  improvement  may  not  be  sufficiently  profit- 
able under  the  circumstances  to  warrant  the  extra  geometrical  and 
other  complexities  involved.   Nevertheless,  the  two  performance 
limits  developed  above  have  considerable  theoretical  interest  and 
value. 

Incidentally,  recall  that  another  simple  but  important  result 
of  the  present  paper  is  its  identification  of  the  key  dimension- 
less  scaling  parameter   X  .   Eqs.  (2.27)  through  (2.29)  show  the 
basic  significance  of  this  quantity. 
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